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Abstract 

Classically, a single weight on an interval of the real line leads to 

moments, orthogonal polynomials and tridiagonal matrices. Appro- 

^ I priately deforming this weight with times t = (ii,t2, ...), leads to the 

Cn ■ standard Toda lattice and r-functions, expressed as Hermitian matrix 

Q . integrals. 

0^ I This paper is concerned with a sequence of t-perturbed weights, 

^^ ' rather than one single weight. This sequence leads to moments, poly- 

f— ^ . nomials and a (fuller) matrix evolving according to the discrete KP- 

hierarchy. The associated r-functions have integral, as well as vertex 

operator representations. Among the examples considered, we men- 

^ ■ tion: nested Calogero-Moser systems, concatenated solitons and m- 

^ ! periodic sequences of weights. The latter lead to 2rn- + 1-band matrices 

l>( I and generalized orthogonal polynomials, also arising in the context of 
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a Riemann-Hilbert problem. We show the Riemann-Hilbert factoriza- 
tion is tantamount to the factorization of the moment matrix into the 
product of a lower- times upper-triangular matrix. 
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Introduction 

The starting point in the standard theory of orthogonal polynomials is a 
single weight p{z)dz on an interval of the real line. The latter leads to mo- 
ments Hij = {z\z^p), depending on i + j only; in turn, moments lead to 
polynomials Pn{z), defined by the determinant (0.2) below and the spectral 
relation zpn = {Lp)n defines tridiagonal semi-infinite matrices L. An impor- 
tant recent development in this ancient theory is that the perturbed weight 

EDO j 

^ 1 '^ p{z)dz leads to t-dependent tridiagonal matrices L{t) satisfying the 

standard Toda lattice equations; the determinants of the principal minors of 
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the moment matrix are r-functions for the Toda lattice and are representable 
as integrals over Hermitean matrices, as developped extensively in |jl|. 

This paper is designed to show the reader how the introduction of an 
infinite family of weights pj{z)dz, rather than a family z^ p{z)dz generated 
by one weight p{z)dz, leads to a theory having many features in common 
with the classic situation above. The weights lead to "moments" fiy, to a 
semi-infinite moment matrix 171^0, to polynomials Pn{z), as in (0.2), and to 
semi-infinite matrices L of type (0.4) below, defined by zpn{z) = {Lp{z))n- 
We mainly deal with: 
(i) t-deformations (t = (ti,t2, ■■■)) 

^,,.,^, . ^'^'p,{z)dz, t=(ti,t2,...) eC°°, 2;gR, j=0,l,2,... 

of the weights; they imply for the matrix L the so-called "discrete KP- 
hierarchy" in t; this hierarchy is fully described in p|, and a large class 
of solutions is explained in section 1. 

Occasionally, shall we deal with 
(a) {t, s)-deformations Q 



Pj 



(t,s;^) = eSr*«-^^F,(-s)p,+,(^), t,seC^, z eK, j = 0,1,2, .... 



e=o 



of the weights pj] they imply for L the 2d- Toda hierarchy, as described in 
| 16| , ^ and summarized in section 2. 

To be specific, given a family of weights pQ{z)dz, pi{z)dz, ... on R, and 
their t-deformations 

Pj{z)dz := pj{t; z)dz = e^i *'' pj{z)dz, 

define the "moments" , with regard to the usual integration in R: 

Piy.= {z\ pj{z)) and pij{t):={z\p]{z)), (0.1) 

and the moment matrix 

m„(t) := (/^ij(t))o<ij<„_i- 

Then the semi-infinite moment matrix m^ satisfies the linear differential 
equations 



dmoo 



dU 



"001 



-'^where the Fi{t) are the elementary Schur polynomials e^i '^ — 'YliiLo Fiit) z^ 
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where A denotes the standard shift matrix. They form a infinite set of 
commuting vector fields. Generically the semi-infinite moment matrix rrioo 
admits a (unique) factorization into upper- and triangular matrices 5*1 and 
S2 respectively, with 5*1 having I's on the diagonal: 

rrioo = Sf ^5*2. 

Consider the vector p{t,z) := {Pn{t,z))n>o of monic polynomials in z, 
dependingp] on t = (ti, ts, •••) e C°°, 



/ Aioo(i) 



Pn{t,z) := Six{z) :-- 



det rrinit) 



det 



P-0,n-l{t) 



f^n-l,o{t) ■■■ f^n-l,n-l{t) 



\ l^no{t) 



/^n,n-l(^) 



1 \ 



,n-l 



2" / 



The eigenvalue problem 



zp{t, z) = L{t)p{t, z) 



(0.2) 



(0.3) 



or, alternatively, the S'l-matrix in the factorization above, gives rise to the 
semi-infinite matrix 



= A + aA° + A^6 + A^^c + .. 



/ ao 1 \ 
bo ai 1 
Co 61 02 1 



v 



(0.4) 



/ 



The polynomials pn{t,z) also give rise to a Grassmannian flag of nested 
infinite-dimensional planes ... D W* D W*_,_]^ D ... , given byQ 

7i:+D W* :=W„e-^i *^^ :=span{p„(t,z),p,+i(t,2;)...}. (0.5) 

We shall also need the associated "Vandermonde" determinants^, 

Ai^)(^) = det(p,_i(2:,.)), <,,<, , ^n{z) = det(4-i)i<,,,<„, (0.6) 

^X{z) ■■= (1, z, 2^ ...) and x*(2) := x{z^^)- 
^'H+ := span {l,z,z^, ...}. 
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and the simple vertex operator, 

..,.,., . . 1 t»^e ^1 ~^; (0.7) 

this is, in disguise, a Darboux transform acting on KP r-functions. We now 

state: 

Theorem 0.1 Given the moments (0.1) and the construction above, the 
semi-infinite matrix L in (O.4) satisfies the discrete KP hierarchy 

BT 

— = [(L")+,L], n=l,2,..., (0.8) 

and has the following r-function representatior^ 

L = td^ag( ^^^^^'--'^^°'A A^^ (0.9) 

e=o \ Tn+2-ern /„>o 

in terms of a sequence of r-functions (tq = 1, ri, r2, ...), which enjoys many 
different representations: 

Tn{t) = det (/if^fc(t))Q<^^^^_-^ (moment representation) 

(0.10) 

1 /■ /■ " 

= —:/•••/ ^niz)A^\z) TT (e^*'^*dzfc) (integral representation) 
n! J JR" ^J--^ ^ ' 

(0.11) 
= det (Proj: e~^i *'^'2;~">V„ ^ 7i+j (fiag representation) 

(0.12) 
= det (Proj : e^i *'^ z"W* — ^ TY+j (dual fiag representation) 

(0.13) 

= ( / (i^; 2;"^"'^p,„_i(2;)X(t, 2;) j 7:„_i(t) (vertex representation) 

(0.14) 

^ where Flip) = Fe{-^, ^^, 5^1 ■••), foi' the elementary Schur polynomials F(. The 
symbol Fi{d)f o g '\s the customary Hhota symbol. 
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where 

"u^ 1 '^ = span{p„(t,2;),p„+i(t,z)...} C 7^+ 

W^e^r*'-' = span||y ^^' _^ " , j = 0, .■., n - 1 1 ® 7^;+ J D 7i;+. 

(0.15) 

The polynomials (0.2) have the following representations 

Pn(i,2;) = JTT TTxx ^^^^ (0.16) 

aet(^/ijj(^tjjo<i,j<n-i 
Ace " 

n\Tn[t) J 7r" ^^^^^ V / 

(0.17) 

and satisfy the eigenvalue problem Lp = zp. 

Notice that formulae (0.16) and (0.17) go in parallel with (0.10) and 
(0.11). Formula (0.17) is a generalization of a formula for classical orthogonal 



polynomials already appearing last century in the work of Heine [11 



We shall apply this theorem to a variety of examples, corresponding to 
sections 4 to 8 (5(x) is the customary delta-function, i.e., J-^6{x)f{x)dx = 

fm ■■ 



Pj{z) := Z-' p{z) tridiagonal matrix L 

Pj+km{z) := z^™'pj{z) 2m + 1-band matrix L 



m 



Pk{z) = S{z — pk+i) — \'l_f_i5{z — Qk+i) concatenated solitons 

Pk{z) = 6'{z — pk+i) + ^k+iS{z — pk+i) nested Calogero-Moser systems 

Pk{z) = {—l)^5^^'{z — p) — 5^'^\z + p) upper-triangular L^. 



The first example leads to the standard Toda lattice and the the classic 
theory of orthogonal polynomials. Since the work of Fokas, Its and Kitaev , 
the Riemann-Hilbert method is a device to obtain asymptotics for orthogonal 
polynomials; for semi-classical asymptotics, see Bleher and Its |0. We show 

Riemann-Hilbert factorization <^==^ factorization rrioo = 5^ ^5*2. (0.18) 
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To be precise, we show the Riemann-Hilbert matrices Yn take on the following 
form (x(^) and x*{z) are as in footnote 2 and hn-i := Tn/T„_i): 



{SMz))n M^2X*W)„ 



YJz) 



Tn (t) Tn (t) 

r„_i(f-[^-i]) ^n-l r„(f+[z-i]) ^-n 

Tn (t) T„ (i) 



(0.19) 

The second example, which is novel and which is developped in section 
5, involves a finite set of weights 

Po{z)dz,...,pm-i{z)dz 

for m >2, which we extend into an infinite ^^m-periodic" sequence 

Po{z)dz, ..., pm-i{z)dz, z"'po{z)dz, ..., z"'pm-i{z)dz, 

z^"'po{z)dz, ..., z'^"'p„,^i{z)dz, ... . 

This sequence leads naturally to generalized orthogonal polynomials Pn{z) 
by the recipe (0.2), which enjoys the following properties: 

i) the polynomials Pn{z) satisfy the orthogonality relations 

{pi{z),pj{z)) = for i > j + 1; 
ii) Applying z"^ to the vector p{z) := {po{z) , pi{z) , ...)leads to a 
2m + 1-band matrix L™. 

„.^ ^^^^ „ ^.^^ .„.. ^ 1 '^ Pk{z) implies L evolves according to the 

discrete KP hierarchy. 
The discrete KP-hierarchy on 2m + 1-band matrices has been studied in 
T7|; see also |lO]. We also formulate here a Riemann-Hilbert problem, which 



should characterize the generalized orthogonal polynomials. 

Another interesting set of examples is provided by picking as weights 
various combinations of (standard) (5-functions, which lead to concatenated 
soliton solutions, Calogero-Moser systems, etc... 

We wish to thank Leonid Dickey for insightful comments and criticism, 
which lead to recognize the importance of the integral in (0.15) beyond its 
formal aspects. L. Dickey has shown in a very interesting recent paper p[ that 
the discrete KP hierarchy is the most natural generalization of the modified 
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KP. We also thank Alexander Its and Pavel Blelier for having explained to us 
the Riemann-Hilbert problem, and for having posed the problem of finding 
the connection with the matrix factorization of the moment matrix. We also 
thank Taka Shiota for a number of interesting conversations. 

1 Vertex operator solutions to the discrete 
KP hierarchy 

In 0, we discussed the discrete KP hierarchy and found a general method 
for generating its solutions, in both, the bi- and semi-infinite situations; this 
paper mainly deals with the semi- infinite case. In [Q and [Q, we gave an 
application of the bi-infinite discrete KP to the g-KP equation. In general, 
the main features are summarized in the following statement, whose proof 
can be found in 0: 

Theorem 1.1 From an arbitrary KP r- function and a sequence of real func- 
tions {...,h'_i{X),i'Q{X),h'i{X), ...), defined on R, one constructs the infinite 
sequence of r-functions: Tq = r and, for n > 0, 

ro(t) = r(t) 

r„(t) = (^y'x(t,A)z/„_i(A)rfA...y'x(t,A)z/o(A)rfA)r(t), n>Q 

r_„(t) = (f X{-t,X)u_n{X)dX ... f X{-t,X)u_i{X)dx)T{t), n > 0. 

(1.1) 
Then the bi-infinite vector 

nt,z)- M'-^'''K ^Tu^^A (1.2) 

V ^nW JneZ 

and bi-infinite matrix 

L = ± d^ag 1^1^,11^1^] A^^^ (1.3) 

i=0 \ Tn+2-eTn J ^^^ 

satisfy the discrete KP-hierarchy equations for n = 1, 2, ....■ 

d^ dL 

— = (L'=)+^ and — = [(L'=)+,L], wtth L^{t,z) = z^!{t,z). {l.A) 

Otk Otk 

8 
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Then Tn{t) is given by the following projection 

Tn{t) = det (^Proj : e'^^ '^'^ z-'^Wn — > 'H+) , (1.5) 

where the Grassmannian flag ... D >V„ D W„+i D ... is given by 

Wn := spand^nit, z), ^„+i(t, z), ...}. (1.6) 

Conversely, a Grassmannian flag ... D >V„, D W„+i D ..., given by (1-6), with 
functions \l/„(t, z) satisfying the asymptotics \[^„(t,2;) = e^*'^' z^ {1 + 0(1 /z)) 
leads to the discrete KP-hierarchy. 



Remark : A semi- infinite discrete KP-hierarchy with T^it) = 1 is equivalent to 
a bi-infinite discrete KP-hierarchy with r_n(t) = Tn{—t) and To{t) = 1; in the 
above theorem, this amounts to setting To{t) = 1 and z/_n(A) := z/„_i(A), n = 
1, 2, .... We extend the semi-infinite fiag Wo = 7i+ D ... D W^ D >V„+i D ..., 
by setting >V_„ = VV*, for n > 0. 



2 Moment matrix factorization and solutions 
to discrete KP and 2d-Toda 

In (0.1), we considered t-deformations of the sequence of weights, with t e 

C°°, 

ipliz),p\iz), ...) , teC^ with p]iz) = e^^'^-'p.iz). 

As announced in the introduction, we consider further deformations of the 
sequence of weights, in the (t, s)-direction, 

oo 

p,(t,s;z)=eSr*«-'^F,(-s)p,+,(z), t, s G C°°, ^ G R, j = 0, 1, 2, .... 

(2.1) 
and the corresponding moment matrix 

mn{t,s) = {pij{t,s))o<i,j<n-i, with pij{t,s) = {z\pj(t,s;z)). (2.2) 

We now state the following proposition (e.g., see 0, ^): 
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Proposition 2.1 The matrix moo{t,s) satisfies the differential equations 

^ = A'm^, ^ = -m^A'\ (2.3) 

dtk dsk 

Factorizing the matrix mao{t, s) into the product of lower- and upper-triangular 
matrices Si and S2, with Si having 1 's along the diagonal: 

m^it,s) = Si\t,s)S2it,s), (2.4) 

the sequence of wave function^ derived from Si and S2, 

^i{t, s; z) = e^^^^^Saiz) mt, s; z) = e-^^^^\Sj)-\*{z), (2.5) 

can be expressed in terms of r-functions Tn{t, s) = det rrin, as follows: 

V Tn{t,S) JneZ 



SiZ 'n 



^2{t,s;z) = — e^i 

V rn[t,s) /„ez 

V rn+l{t,s) ;„g2 

nifs;z) = f^MlMe-Er --.-") , 



(2.6) 



with \l/j(t, s) satisfying the following differential equation^ 

-^ = m)+^i , -^ = {^)-'^^ with Li = SiAS^\ L2 = S2A-'S,\ 
The T-functions satisfy bilinear identities, for all n,m >0, 

T^{t-[z-'],s)T,^+i{t'+[z-'],s')e^T(^'-^>'z--^-^dz 



r„+i(t, s - [z])TUt', s' + [z])e^T(^^-<>-'z^-"'-'dz, (2.7) 

and therefore the KP hierarchy in each of the variables t and s. 



z=0 



^ ^i{z) ■■= Y.Uz' and ^2(2) ■.^Y.SiZ *; also x{z) ■= {^, z,z'^, ...) and x*{z) -^ x{z ^)- 
"^ A-i- and A^ denote the upper-triangular and strictly lower-triangular part of the ma- 
trix A, respectively. 

10 
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The following corollary can be found in |^: 

Corollary 2.2 2d- Toda t -functions satisfy the following (Fay-like) identities 
for arbitrary z,u,v G C. 

r„(t - [Z''^],S + [v] - [u])Tn(t, S) - Tn{t, S + [v] - [u])Tnit - [z^^], s) 

= -^Tn+iit, s - M)r„_i(t - {z^\ s + H). (2.8) 

Introduce now the residue pairing about 2; = oo, between / = X]i>o '^i^* ^ 
H+ and g = Ejez bjZ'^'^ e H: 



{f,g)oo=i f{z)g{z)^ = J2aA, (2.9) 



where the integral is taken over a small circle about z = 00. 

But setting s = s' and m < n — 1, the right hand integrand of (2.7) is 
holomorphic and so the right hand side of (2.7) vanishes. Of course, freezing 
s = s' yields the discrete KP-hierarchy; see [Q. Therefore 

/ Tn{t-[z-'],s)Tm+iit'+[z-^],s)e^^^'^-''^^'\'""'^'dz = ioi n> m + 1, 

(2.10) 
and so for n > m -\- 1, 

I ^T.TUz^.n rn{t-[z-'ls) ^^y^,,^.^_^_, T^^,it'+[z~'],s) ^^ ^ ^ 

Jz=oo Tn{t,S) Tm{t, S) 

(2.11) 
Defining the linear space W* as the space of functions perpendicular to >V„ 
for the residue pairing (2.9), we thus have for fixed s, by virtue of (1.6), (2.6) 
and then (2.11), 

W* = span{;.^-^^^^^i^4li^, j>n} = e-^'^^'Wn 

Tj{t,s) 

W:* = span{z-^'^^^^^i^4^' J < n} = e^'^^'W:. (2.12) 

r,_i(t, s) 

It also shows that r„(t, s) can be obtained from those spaces in two different 
ways (for fixed s): 



r„(t,s) = det (Proj : e-S*'"";2-"W„ — > H^ 

= det (Proj : e^'^'' z'^W^ — > 71+) , (2.13) 



11 
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where the multiphcation by z~" and z" makes the corresponding hnear spaces 
have "genus zero", in accordance with the terminology of Segal- Wilson |jT4|. 

As a special case (Hdnkel matrices) , consider the sequence of weights 

Pj{z)dz = z^ p{z)dz. (2.14) 

Then the (t, s)-deformations take on the following form: 

p,(t, s; z) = e^*'^' Y. Fe{-s)z'+^p{z) = e^^'^-'^>' z^ p{z), 

£>0 

thus depending on t — s only. Therefore Pij{t, s) depends only on t — s and 
i—j [rrioo is a Hankel matrix) and so r„(t, s) depends only on t — s. Therefore, 
in this case we may replace t — s hj t. 

In this case, the matrix rrioo is symmetric, which simplifies the factoriza- 
tion (2.4) above. Indeed: 

moo(t) = S^^S2 = S^'hS^'^ = S-\t)S^-\t), (2.15) 

upon setting 

S = h-'/^S^ = h'/^S,'^. (2.16) 

3 Weights, flags and dual flags 

The purpose of this section is to prove Theorem 0.1. The point is to derive 
the r-functions from the Grassmannian flag (1.6). Unfortunately, the matrix 
associated with the projection (1.5) is infinite; therefore taking its determi- 
nant would be non-trivial, although possible. However, it turns out to be 
inflnitely easier to consider the dual flag, which leads to a finite projection 
matrix, whose determinant is the same r-function. 

To carry out this program, we equip the space H := span{z*, z G Z} with 
two inner products: the usual one 



{f,9)= / f{^)9{^)dz, (3.1) 

and remember the residue pairing about z = oo, between / = I]j>o '^j-^* ^ '^^ 
and g = J2j£Z bjZ~^~^ E H: 

^^' 3)o^ = i_^ f^'^3{z)^^ = E «^&^, (3-2) 



i>0 



12 
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where the integral is taken over a small circle about z = oo. The two pairings, 
which will be instrumental in linking the flag to the dual flag, are related as 
follows: 



Lemma 3.1 



{f,9) = (f,[ ^^du) . (3.3) 

\ J-R Z — U / 



Proof : Expanding the integral above into an asymptotic series, which we 
take as its deflnition, 



^^du 



R Z U ^ ^±x. j^Q 



-Y.^'' I 9W)u'du, (34) 



/, / du) . = /^flj^;*, -^2; ^ I g{u)u^du, 



we check that for holomorphic functions / in C, 

00 \j>0 ^ j>0 

^Oj / g{u)u^du 
g{u) y^ aju'du 

i>0 

{f\g). (3.5) 



Remark : The series (3.4) only converges outside the support oi g{u). So, 
in general, the series (3.4) diverges, even for large z. In speciflc examples, 
this integral will have a precise meaning; see sections 4 and 5. 

To the family of functions po(-2)) Pi (-2)) ••• on R, andp*(z) := 6^**=^ Pj{z), 
we associate the flag of spaces Wq = 7i+ D ... D W„ D W„+i D ..., deflned 
by 



Wn ■■= (span{pi,pi,...,p„_i}^ 

= {f en+ such that (/, Pi) = 0, 0<i<n- 1} (3.6) 



13 
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with respect to the inner product (3.1). So, throughout we shall be playing 
with the following two representations of the moments: 



/ii,- = {z' , p]{z)) = iz' 



p]{u)du^ 



R Z — U 



(3.7) 



With the moments Pijit) := {z\p'j), we associate the monic polynomials 
Pk(t,z) in z of degree k, introduced in (0.2). As usual, set 

W*=e-S*»^'Wn and its dual W** = e^^^^'W^. 
As we showed in (2.12), for the residue pairing we have: 

The integral representation (3.9) below of the dual flag already appears in 






the work of Mulase |jT3|, for the case pj{z) = z^ p{z). 



Proposition 3.2 The flag n+ D Wi D W2 D ..., defined by (3.6) att = 0, 
evolves into 

W* = (span{p*,p*,...,p*„„J)^ = span{p^{t,z),pn+i{t,z),...} C n+, (3.8) 

and the dual flag 7Y+ C W^ C W2 C ..., evolves into 



W:'=span(U 



pUu)du 



z — u 



j=0,...,n-l ©T^H 



(3.9) 



Proof : Indeed to show (3.8), it suffices to check the following, for k > j + 1 

k 

and the polynomials pk(t,z) = ^ f^^ y^aki{t) z\ defined in (0.2): 



j=0 



{pk{t,z),p') 



akkif) 

1 
akk{t) .^0 



^afci(t)(z\p5) 

i=0 
k 



( poo(i) ••• l^O,k-l{t) 



akk{t) 



det 



\ pko{t) ••• Pk,k-l{t) 



f^Ojit) \ 



0. 



f^kj{t) J 



(3.10) 
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To prove the dual statement (3.9), one checks for A; > j + 1 

\ JR Z — U 



oo 



using Lemma 3.1, and, of course, 

[pk{t,z),z^) =0, for all A;,£> 0. 



Remember from (2.13), the r-functions r„(t) can be computed in two 
different ways: 



r„(t) = det(Proj :e-S*>"'z-"W„ — > H^ 

= det (Proj : e^'^'^ z^'W;^ — > n+) . (3.11) 



We shall need the following lemma concerning Vandermonde-like determi- 
nants, extending a lemma mentioned in [|T^]: 

Lemma 3.3 

^ det (u^,a(k)Vk,a{k))^^^^^^ = det (w^,fc)i<^,fc<„det {ve,k)i<e,k<n ■ (3-12) 



Proof of theorem 0.1: Since z^W* D -z"7i+, the matrix of the projection (3.9) 
onto TC+, involving W*, reduces to a finite matrix, whereas the projection 
involving Wn would involve an infinite matrix! This is the point of using W* 
rather then W„. Therefore the matrix of the projection 

Proj : eS*^^'2"y\;** — > n+ 
is obtained by putting all coefficients of 

eE*^^";^" f Pi(^)du ^^^ (0<j<n-l) and eS*^"'z"+^' for (0 < j < oo) 
J z — u 

in the jth and n + jth columns respectively, starting on top with 2;°, 2;^, .... 
Since for any power series 

z^-coef of / = / z-^-^f{z)^. = {z~^~\ f{z))^, 

Jz=oo ZiTt 
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we have 



where 



r„(t) = det 



A 



B 



C 



det A det C = det A, 



(3.13) 



C = (coef,n+.z"+^eS*'=^ 



0<ij'<co 



/I ... \ 
Fi 1 ... 
F2 Fi 1 ... 



V : 



•• / 



and 



A 



, V J-R. Z — U 



^n-i-lg. 



Y^t^z^ f Pj{u)du ' 



R Z — U 



0<i,j<n—l 



00 / 0<ij'<n— 1 



U 



i-i-1 X]*fe"* 



Pj[U, 



0<«J<n-l 



= (/i„-,-ij(t))o<.^^.<„„,, (3.14) 

which provides the A-matrix in (3.13), thus estabhshing (0.10). Hence, 

r„(t) = det (^^fc(t))o<^;.<„_^ 
= det I / z^pk{t,z)dz] 



0<e,k<n-l 



R 



det ( / z^(^-.pk--i{t, Za(k))dz^tk)] for a fixed permutation (T 

^ '" ^ ' J i<e,k<n 



= -/•••/ det(4-^)i<fc,^<„det (p£-i(2fe))i<fe,^<„ n (e^*''^c?Zfc 
using Lemma 3.3; this estabhshes (0.11). Furthermore, we have, continuing 
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the identities above, that^ 

JR 

/ A„_i(2;i, ..., Zn-i) n (p£-i(-2^)e^*"^^rf2;^ 

= (^^ rfz^"-Vn-iWX(t,z))r„_i(t), 

proving (0.14). Therefore, the sequence r„(t) satisfies (1.1) in Theorem 1.1 
with Un{z) = z'^pn{z) and ro(t) = 1. The r-functions lead to the expression 
(1.3) for L and to the expression (1.2) for \1/, which both satisfy the discrete 
KP hierarchy, according to Theorem 1.1. Notice that, from (0.15), (1.5) and 
(1.2), we have 

.„v^,~. ^ i*'^^„(t,^); (3.15) 

therefore L defined by r-functions (1.3) agrees with the semi-infinite L, de- 
fined by the semi-infinite polynomial relations zp{t,z) = Lp{t,z), yielding 
(0.8) and (0.9) for this L. 

Finally, using (0.10) and (0.11), the wave function (1.2) equals, 

Tn{t) 

= jTT TTT^ ^^^ , usmg (0.10), (0.1), footnote 8, 

aet [fiij [t))o<ij<n-i 
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establishing (0.16) and (0.17). ■ 

In the subsequent sections, it is shown that many integrable solutions, 
when linked together, are nothing but special instances of the situation de- 
scribed in section 2; we mention matrix integrals, 2m + 1-band matrices, 
soliton formulas, the Calogero-Moser system and others in subsequent sec- 
tions. 

4 Toda lattice, matrix integrals and Riemann- 
Hilbert for orthogonal polynomials 

Setting 

Pj{u)du := u' p{u)du, p {u) := p{u)e^i *=" , (4.1) 

define the moment matrix 

moo(t) := (/iiiW)o<^j<oo with /i,,(t) = J z'+^e^^''^'" p{z)dz, (4.2) 

and the corresponding t-dependent monic orthogonal polynomials Pn{t, z) in 
z. Note that moo is a Hankel matrix and is therefore symmetric. From the 
form of the moments, the matrix m,oo{t) satisfies the the following differential 
equations 

^ = A^moo. (4.3) 

Refering to the special case of Hankel matrices, discussed at the end of section 
2, we consider the factorization of the symmetric matrix m,oo{t) into the 
product of a lower- and upper-triangular matrix 5*1 and 5*2, with I's along 
the diagonal of 5*1 and /I's along the diagonal of 5*2 : 

m^oit) = S^^S2 = S{^hS{^^ = S-\t)S^-\t), with S = h-^'^Si = h^^^S^^^. 

(4.4) 

Theorem 4.1 Then S{t) and the tridiagonal matrix L{t) = S{t)AS^^(t) 
satisfy the standard Toda Lattice equation^: 

f)C 1 AT 1 



^with regard to the splitting of ^ G gloo into a lower-triangular Aj, and skew-symmetric 
matrices Ask- 
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The flag and dual flag of (0.15) take on the following form 

Wl^ = span{pnit,z),pn+iit,z),...} 

= span{(5(t)x(^))„,(5(t)x(^))„+i,...} 

= span < z , < n < oo > 

I ^n(t) - J 

^^,*-^ i fpj(t,u)p\u)du . ^ -,1-1, 

W** = span<^ / ^n > Ji^K J — ,j = 0,...,n-l>®n+ 



span 



{z~\{Sm^{t)) x*{z))j,J = 0, ..., n -l}®n, 



= span|^-^-^ ^^+^^^|J^''J\ j = 0,...,n-l|®7^+, (4.6) 

with the T-functions having the following representation, derived from (0.10) 
up to (0.13), 



Tn{t) = det ( f z'+^ p\z)dz) 



0<i,J<n-l 
1 „ „ n 



= / e'^^(^W+i:r*^Af')rfM, setting p{z) = e""^^^ 

= det [Proj : e" ^^ '^'^ z'^'Wn ^ :H;+) 

= (^l^dzpiz) z^^^-'^Xit,z)yn^,it), (4.7) 

anc? it/ie orthogonal polynomials, having the form 

Tn(t) 

/^^ det(^/ - M)e^'-(^(^)+i:r*»^')rfM 

where dM = Af^{z) dzi...dzn dU is Haar measure on the set of Hermitean 
matrices Tin- 

Before stating the corollary, some explanation is needed. The integral in 
the matrix below is taken over the R with a small upper semi-circle about 
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z, when Qz > and over R, with a small lower semi-circle about z, when 

Qz < . Moreover Yn±{z) = lim ^/_^ Yn{z'). 

±3;z'>0 

Corollary 4.2 /n mew of the factorization moo(t) = Si^S2 of the moment 
matrix m^{t) and setting hn = Tn+i{t) / Tn{t) , we have the following identity 
of matrices: 



Yn{z) 



h-\pn^,{t, Z) h-\ J^ ?^^^p\u)du , 

Un-1 (5lX(^))„-l /^n-lM^2X*W)„_l, 



rn(t-[z-^]) 



z" 



Tn + l{t+[z ^]) ^-n-1 



T-„(t) 



Z 



Tn-i(t-[z~'-]) ^n-1 r„it+lz~'']) ^-n 

T„it) "^ T„(t) ^ 



(4. 



The matrix Yn satisfies the Riemann-Hilbert problem of Fokas, Its and Kitaev 

&■■ 

1. Y{z) holomorphicf^ on C+ and C_. 



2. Y(z) = Yjz) 



1 27rip\z) 
1 



3. Yiz) = iI + Oiz~')^'^^ ^l^^whenz 



CXD. 



(4.9) 



Note the first column ofY{z) relates to the Grassmannian >V„ and the lower- 
triangular matrix Si, whereas the second column to the dual W* and the 
upper-triangular matrix S2 in the decomposition of m^ = S'f^S'2. 

Proof of Theorem 4.1 The vertex representation (4.7) of r„(t) shows that 
the r-vector r(t) = (r„(t))„>o is a solution of the discrete KP equation (1.4). 
But more is true: L = SAS~^ is tridiagonal; so, S and L satisfy the standard 
Toda lattice (4.5). Some of the arguments are contained in [|I|. 



10 



C+ and C_ denote the Siegel upper- and lower half plane. 



20 



Adler-van Moerbeke:Polyn/RH-problems Aug 24, 1998 §4, p. 21 

Notice that the Borel decomposition (4.4) is tantamount to finding the 
orthogonal polynomials pn{t,z) with respect to the inner-product {z'^,z^) = 
fiij, to be precise: 

moo = S-'S^-' ^^ Sm^S^ = I ^^ (/^^'^'p., /^7'^'p.) = ^^3■ (4-10) 

— 1/2 

It follows that the coefficients of the orthonormal polynomials h^ Pi are 
given by the ith row of the matrix S{t) and so 

S^{t) = h^'^S{t) = (p.,(t))o<ij<oo, where p„(t) = ^ Pnj{t)z^ . (4.11) 

0<j<n 

(i) So, the monic polynomials Pn{t,z) of (0.17) have the following form: 

hl/HS{t)x{z))^ = {S^{t)x{z))^ 

= Pnit,z) 

^n rn{t-[z-^]) 

Tnit) 
1 r r " 

= ;^M^ / - L ^'^"^ n ((- - u,)e^^'-hiu,)du,) , 

(4.12) 

leading to the formula in the statement of Theorem 4.1. The p„'s are the 
standard monic orthogonal polynomials with regard to the weight p^{u) = 
p(M)e5I)**"\ 

(a) But, we now prove 

h'J' {S^-\t)x*{z)) = iS,{t)x*{z))^ = z f P-^^^^l^^^^P^du 

^ ^ ' '' J z — u 

^_„r„+i(t+[2;-i]) 



Indeed, we compute, on the one hand. 



. (4.13) 



j>0 j>0 






Zl^ ^ZlP«^W/^<?i' using (4.11) 
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^>o i>o 

Pn{t,u)p\u) 



R Z — M 



On the other hand, as we have seen in the special case following (2.14), the 
2d-Toda r-function r(t', s') depends on t = t' — s' only, enabling us to write 
(here ip stands for \1' without the exponential). 



hl/'j:{s^-\t))^,--' = {h^^'S--\t)xiz--))^ 



n-i 
j>0 



[S2{t',s')x{z-'))^ , using (4.4), 

ip2,n(t',s';z~^) , using (2.5), 
r.+i(t',.'-[^-i])^.„ 



rn{t',s') 

Tn+l{t+[z-']) 

z , 



using (2.6), 



Tnit) 

from which (4.12) follows, upon using Sm^ = S~^~^; see (4.9). 

Theorem 4.1 is established by remembering Proposition 3.2 and using 
(3.8) and (3.9); i.e., 

W* = Span{prr{t,z),Prr+l{t,z),...} 

,*t _ __ [ /■ u^p\u)du 



Wr = spanjy ^^^ ,j = 0,...,n-l|®7^+ 

= span|/ ^^-(^-;y^"^'" ,,=0,...,n-l}e7^,, 

together with (4.12). ■ 

Proof of Corollary 4.2: Following the arguments of Bleher and Its 0, 
the first matrix in (4.8) has the desired properties taking into account the 
following integrals: 

— - hm / — p [u)du = Pn{t,z)p [z) + — - lim / — p {u)du. 

'32'<0 '32'>0 
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The formulas (4.11) and (4.12) lead to the desired result. ■ 

Remark : From the fact that det Yn- = det Yn+, it follows that det Y{z) is 
holomorphic in C and since det Y{z) = l + 0{z~^), it follows from Liouville's 
theorem that detY{z) = 1, i.e., 



detYn = /i„„i p„(t, z) / p {u)du-pn-i{t,z) p{u)du] 

\ Jn. z — u Ju z ~ u I 

= 1. (4.14) 

This is not surprising, in view of the fact that the first expression for det Y^ 
is nothing but the Wronskian of the the two fundamental solutions of the 
second order difference equation; see Akhiezer 0. The second expression, 
involving r-functions follows also from Corollary 2.2, by setting u = z~^ and 
V ^ and by using the fact that, for the standard Toda lattice, we have 
r(t, s) = r{t — s). 

5 Periodic sequences of weights, 2m + 1-band 
matrices and Riemann-Hilbert problems 

The results of section 4 about tridiagonal matrices will be extended in this 
section to 2m + 1-band matrices. As usual, we set 



J J ^ ^\\ „.;j-T, J ^ ^\ _ oX]*fe^* 



p,j{t) = {z\ p]{z)) , with p]{z) 



Pj{z). 



(5.1) 



In proving and stating the results below, we shall also consider the s- deformations, 
as in (2.1). Here we consider m-periodic sequences of weights po, pi, ..., de- 
fined by 

Pj+km{^) = ^'>,(^), for all J = 0, 1, 2, .... (5.2) 



Theorem 5.1 For the weights (5.2), the polynomials 

I IIoo{t) ... Po,n-l{t) 



Pn{t,z) 



det (/i£,fc(t)) 



■det 



0<e,k<n-l 



fJ-n.-lfiit) ■■■ Pn-l,n-l{t) 



\ l^noit) 



fJ'n,n—l\J^ ) 



1 \ 



^n— 1 



2" / 
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= z"'— -^ — , where Tn{t,0) = det rrinit) 

rn[t,0) 

det {zfiij{t) - /ii+ij(t))o<^^j<„„i 

det [fiij (t) )o<j,j<n-i 
Iff " 

= ^^/•••/ Kiz)Ai^\z)U{^^^'''Hz-Zk)dzj:), (5.3) 

lead to matrices L, defined by zp = Lp, 

(i) which evolve according to the discrete KP hierarchy 
(a) such that L™ is a 2m + 1-band matrix . 

(Hi) the polynomials Pn{z) satisfy the generalized orthogonality relations 
{Pi{z), pj{z)) =0 fori>j + 1. 



Remark : It is interesting to point out that the condition (5.2) is equivalent 
to a seemingly weaker one: 



z'^pj e spanjpo, ■■■,pm+j}, for all j = 0, 1,2, ..., 



(5.4) 



where Pm+j must appear in the span. Indeed, the p„'s only depend on the 
moments pij by means of the determinantal formulae (5.3), which allow for 
column operations. 

Corollary 5.2 The following 2x2 matrices are all equal 
I 



Yn[z) 



Pn{t,z) 



p..(M) fl^.-.^_^^^,^ \ 



\k=l 



K'-^Pn^lit, Z) h-U i ^^^^^ it Z-'pU{u)] du 



\k=l 



iSMz))r 



\{s,r{z)\ 



( r„(t-[2;-i],0)^„ r„+i(t,-[2:'i]) ^_„_,\ 

rn(t,0) 



r„(t,0) 



V r„(t,0) 



r„(t,0) 



/ 



(5.5) 
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they solve the following Riemann-Hilbert problem: 
1. Yn{z) holomorphic on C+ and C . 



2. Yn-{z) = Y,,+ {z)(y 



1 



5. y„(z)(^^Q ^nj — >I,asz^oo. 

4- Yniz) L ^^_i j finite , as z^O. 

The polynomials Pn are such that z"^Pn satisfies 2m + 1-step relations. 
Note 

"" r„(t,0) 

and the first column of ¥„ is related to the Grassmannian plane W„ and the 
second column to a plan^ related to ^2it^ s, z). 

Remark : In the matrix (5.5), r„(t — [z~^],0) is given in formula (5.3), 
whereas by (0.10) and (5.11), 

r„«. -[.-']) = det((«'.eEr..-'g£i;^\) . (6.6) 

V\ fe=0 ^ " I ^ 0<i,j<n~l 

Also note the right hand column of (5.5) behaves as 1/2"+^; this follows from 
the r-function representation, but also from the "generalized orthogonality", 
mentioned in (iii) (Theorem 5.1). 

Proving Theorem 5.1 requires the following lemma: 

Lemma 5.3 Fix m > 1; the polynomials pn{t,z), defined in (0.2), satisfy 
2m + 1-step recursion relations, i.e., 

z'^plt, z) = L"^p{t, z) with L™ a 2m + 1-band matrix, 

if and only if every pj, j = 0, 1, ... satisfies the following requirement: 

For every i = 0,1, ..., j = 0, 1, ... there exist constants Cr, r = 0, ...,m + j + 
depending on j and i such that 

{u'^pj - E™^^^^ CrPr, M^) = /or < i < m + j + ^ + 1 



^^ Unlike the orthogonal polynomial case, the second column does not contain elements 
of the dual Grassmannian W* . 
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Proof : Note the following equivalences: 

z"^Pn{t,z)= ^ A(t)nrPr(t,z) for souie matrix 74(t) 

n—m<r<n+m 



z"'Pn{t,Z) G y^ma.x{n-m,0) 



^ ' "^n ^ ' '^max(ra— m,,0) 



for all n > 

for all n > 0, because 
of the inclusion 



2™>V„ C W, 



max(ra— m,0)- 



Since 



W„ = (span{po,Pi,---,Pn-i})"^ = span{pn{z),pn+i{z),...} 
the latter is equivalent to 

= {u"^pn{u),pj{u)) for all < j < n - m - 1, n > 

= {pniu),u"'pj{u)) 

1 



(^nny^) {—Q 



J20'ni{u\u^Pjiu)) 



(^nny^j 



det 



\ pnO ■■■ Prijn—l 



fJ"mj \ 



Pm+nJ / 



where we have used the fact that Pn(t, z) = ^,^. X) o,ni{t)z^ is represented by 
(0.2). The vanishing of the determinant above is equivalent to the statement 
that the last column depends on prior columns; namely there exist Cq, ..., c„_i 
depending on m, n, j such that 



n— 1 



= Pm+i,j — X! '^rf^ir foTO<i<n, 0<j<n — 171 — 1 
r=0 

71-1 

r=0 
n—1 

= {u\u"^Pj -^^CrPr) 
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j+m+e 

= {u\ u"^pj - ^ CrPr) for < 2 < m + j + £ + 1, 

r=0 

where d, was defined such that j + £ = n — m — 1. ■ 

Proof of Theorem 5.1 : The fact that L™ is a 2m + 1-band matrix follows 
at once from Lemma 5.3. That the matrix L evolves according to the dis- 
crete KP-hierarchy follows straightforwardly from the general statement in 
Theorem 0.1. ■ 

Proof of Corollary 5.2: For the sake of this proof, we shall be using the 
(t, s)-deformations of the weights pj and the corresponding matrix moo{t, s) 
of (t, s)-dependent moments 

oo 

fiij{t, s) = {z\pj{t, s; z)) , with pj{t, s; z) = e^i **^' ^F^(-s)pj+^(2;). 

1=0 

(5.7) 
In section 2, it was mentioned that rrioo satisfies the differential equations 

^ = A'm^, ^-^ = -m^K'\ (5.8) 

otk dsk 

Factorizing the matrix 

m^{t,s) = S^\t,s)S2{t,s) (5.9) 

into the product of lower- and upper-triangular matrices 5*1 and 5*2 then leads 
to the 2d Toda lattice. 

For later use, we also compute pfe(t, — [z~^]; -u), making specific use of the 
periodicity of the sequence of weights pj^rn{u) = u^Pjiu) and the identity^^ 

F^{\z'^])=z-^. (5.10) 

We find: 

oo 

pu{t,-[z-\u) = ei:r*^«'^F,(ri])p,+,(«) 

j=0 



^^obtained, by expanding the following expression in elementary Schur polynomials, by 
setting t ~ Q and by comparing the powers of y: 



^,«F„(t+[.-])=en*-^y=eE*--(i-^)-^.5:,«^„(t)f:m\ 



n>0 



27 



Adler-van Moerbeke:Polyn/RH-problems Aug 24, 1998 §5, p. 28 



e2.i *'« Y, ^-^^^, using again (), 



i=o 

'm— 1 



2m- 1 



7=0 ^ i=m ^ 



m— 1 

j=0 
m— 1 

j=0 






2m 



+ ... 






Pj+fciMj 



'M^m^ 



eEr*^«^^ 



m— 1 ^m-j „ /„,\ 



3=0 



Z^ -M™ 



(5.11) 



From (5.9), we have Sirrioo = 5*2 and hitting x*{^) with this matrix, we 
compute, on the one hand. 



i>o 



l"^ooj„j- ^ 



s=0 



J2^ ^J2Pne{t,s)fii 



j>0 i>0 



s=0 



j>0 £>0 ^ 

^^r--'' /" p„(t,s;M)eS"*'"'^F£(-s)pj+£(M)rf^ 
/ p„(t, 0; M)eEr *''-' ^ F,([^-i])p,(iz)dM 



s=0 



R 



Pn{t,0;u)po{x,-[z ^];u)du 



/^P„(t,0;.)eEr.,.'X:£_i^A, (5.12) 



using (5.11) in the last identity. On the other hand, we have 



E(^2(t,s)) 



nj 



s=0 



-1- 



S2{t,s)x{z 



^2,n(t,0;2;-i) 
rn(t,0) 



s=0 



(5.13) 
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The right hand sides of (5.12) and (5.13) coincide (using Sirrioo = S2)', 
so, we find the following identity, together with the desired asymptotics for 

z ^ 00: 

(5.14) 
leading to condition 3. The jump condition 2. follows from the following 



1 ,._ f Pn{t,u 

2TTi 









^^ ' V ^ '^■'^ ' 2711 .'^z J z'^-u^'x^^ I'jy > \ 



(5.15) 
The asymptotics 3. follows from the r-function representation of the first 
integral. 

The formula concerning det Y^ follows from setting m = and v = z~^ in 
identity (2.8). ■ 



6 Soliton formula 

For future use, we define the vertex operator: 

X{t- A, /i) = ^^eS"(^'-'^')*'=e^"^''"'"^"'^^^. (6.1) 

A — 11 

Theorem 6.1 Given points p^, g^ and X^, k = 1, ..., and the weights 

Pk = 6iz - pk+i) - \l+iS{z - Qk+i), /c = 0, 1, ... , 
the T-functions 

rJt) = detfe~'eSr=i*'=P?-A,2gr'eSr=i*'=9n 

V / l<i,i<n 

= {pl-'X{t,Pn) - \lq:-'X{t, qn)) r„_i(t) 
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fc=l V ^« Pi J l<i,j<n 



It IL 

fc=l 1 

/orm a r-vector of the discrete KP hierarchy, for appropriately chosen func- 
tions ai and Cn of p,q,X. The matrix L, constructed by (0.9) from the t's 
above, satisfies 

L{t)p{t,z) = zp{t,z), 

with polynomial eigenvectors (in z): 

det (5iAz-pi) - ^tiiz^SileT.ZM'i'-p')) 
detf(5,;. --r^^eSfe=i*^('?»-*'^) 

Then 



^-^ I^-PJ " I l<ij<n 



>V„ = (span{po,---,/On-i})"^ 

= {/g7^+, suchthat f{p,) = \\f{q.:),% = \,...,n}{^-'^) 
and 

W: = span I ^ ^,i = l,..,nl©7^+. (6.3) 

\z-pi Z-Qi J 

Proof : Consider the space 

n+/z''n+ = span{l,...,2"-i} 

= spaniel,..., w„}, (6.4) 

where the polynomials 

Vk{z) = J] (z-pj) 

l<j<n 

of degree n — 1 form an alternative basis; the determinant of the transfor- 
mation between the two bases being the Vandermonde determinant A„(p) = 
Ili<j<k<niPk-Pj)- Define 

n 

Piz):=l[iz-p,) (6.5) 
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and 

Vk{Pk) ^n[p) 



With this notation 



n 
fc=l 



Vk{Pi) = SikVkiPk) = SikP'iPk), Vk{qi) = —■ (6.6) 

Qi -Pk 

Using the second hne of (3.14) and the formula for pi, one computes 

Tn{t) = det(/ijj)o<ij<ra-l 

= ^'4f '"""««'"" (A - A) Is) 

\J z=oo \Z — P{ Z — C/j / ZTil, J 

0<j<n-l 

using (6.4). The second identity above leads to the first formula in the state- 
ment above about r, whereas (0.14) is responsible for the second formula. 
The last formula on the right hand sidejust above, leads to 

T-n(t) = -^-T det {vj{v^)g{J)i) - Vj{qi)g{qi)\'^i)^^,,^^ 
det diag{g{pi))i<i<n 



A{p) 



S^MP^)9iP^) - ^i (^*)^(9*)^*^) i<ij<„ ^^^g [g{p^) \ ^^^^^ 
I Q ( T) 

det ( SijVi{pi)g{pi) - \}vj{qi)g{qi)-j-!^ \ , using (6.4) 



Cnit) det 5ij - \ 



,2vMi)g{qi) 



V ^M) 9{Pj) J i<ij<n 

c„(t)detU.,-— — -— 

\ yi Pj yypjj / i<i,j<n 

Cnit) det {6ij - aiX{t, qi,pj)l)^^..^^ 

n 

1=1 
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using in the last equality the vanishing of the square of the vertex operator. 

The formula for Pn{t, z) is derived from the third expression for r(t), using 

the standard representation (2.6) for the wave vector \l/(t, z). ■ 

Remark : When g^ = — p,, the formula for the KdV r-function reads: 

(n \^n 



dd "'^i 



Cnit) det [5,, + ^^e~Sodd*^(p.'+P.') ] 

V Pi + Pj J l<i,j<n 



Note that Segal and Wilson have used, in [Q, the infinite matrix represen- 
tation of the projection of (6.2), rather than (6.3), in order to compute KdV 
solitons. 



7 Calogero-Moser system 

Theorem 7.1 Given points pk, Xk (k = 1, 2, ...), the weights 

pk = 5'{z-pk+i) + \k+i5{z-pk+i), k = 0, l,...,n- 1 
determine a sequence of r-functions for the discrete KP equatior^, 

Tnit) = -, - \{e^'^'''^^n{z)^^P\z)dz^...dZn 

= etrTT ^^y- det (-X + f] khY^A , (7.i; 



with appropriate matrices X and Y , functions of pk and Xk's, satisfying the 
commutation relatior^ [X, Y] = le, and having the form 

X = diag(xi,...,x„) and Y = i — — ^ + diag(^i, ...,^„). (7.2) 



"where le = (1 - %)i<ij<r 
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The matrix L, constructed by (0.9) from the t's above, satisfies 

L{t)p{t,z) = zp{t,z), 
with eigenvectors, polynomial in z, 

Pnil z) = det i zl -Y - I xI + f^ktkY''-^-x\ j. (7.3) 

The Grassmannian flag corresponding to this construction is given by 
Wn = {fen+, such that f'ipi) = Xi f{p,), l<t<n}, 
W: ^ {^-^,.^l,...,.)e7^.. (7.4) 

Proof : As before, we introduce the basis Vk{z) of H^ /z^H~^; note 

U^ i = ljVfc,i i=l i^k Pi 

and the matrices 

^ = -diagfi:^ A.) -(i^) 

\ l<a<n Pi Pa j \Pl Pj / l<ij<n 

Y = diag(pi,...,pn) (7.6) 

with commutation relation 

[X, Y] = h, /e = (1 - 5^j)i<ij<n. (7.7) 

Then, by (3.7) and the choice[^of pj, 

r„(t) 

= det(/ijj)o<jj<n-i 

1 A,- \ dz\ 



det ( / z^'^g{z) 

\ J z=oo 



\{z-pi)'^ Z - Pi J 27Ti J ^^..^^ 



^^c„ :— A (' )' ^^ ^^^ expressions below. 
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1 



det f Vj{z)g{z) 



A,- \ dz 



{z-piY z-pij 2m J ^^. 



i,j<n 



A.(p) '"^ I ("'"'' 



-\iVj{pi)g{pi 

^-Pi / l<i,j<n 



i^ det ( gip,) ± ^^'^^-^ ^"^^-^ + v.ipMP^) (t ktupt' - A. 

An(p) i f^ Pi -Pa Vl 



/ 



fc=l 



1 - <5,;, 



nei:r*«Pldet ^^^ + 5,J^ 



v^ 



Pi - Vj 

oc 

etr Er *>^' det ( -X + ^ HfcF^-M , 



-| oo 

2 — + ^ ktkP^;-^ - \ 



l<i,j<n 



l<i,j<n 



yielding the formula for r„(t); According to theorem 0.1, the pn{t,z) are 
polynomials, which we now compute: 



Pn{t,z) = Z 



nrn {i-[z-']) 
Tnit) 



oo / y' 

det ( -X + y fctfcF^-i - z-^y 



fc-lN 



1 



2; 



det [-X + Y, ktkY''-^ 

/ _ oo _ _ 

det (-X + J2 kikY''-^ - z'^ (l - z'^Y 



det{zl - Y) 



~\ -1 



det l-X + y ktkY 



fc-i 



det(z/ -Y) det / - hrJ + ^ ktkY''-^ - x] (z-Y) 



-1 



(7i 



yielding (7.3), but also an expression for the wave functions ^„(t,z), upon 
multiplying by an exponential. The formulae for Wn and W* follow from 
(3.8) and (3.9) and the choice of pk- 
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In order to connect with the form of the matrices announced in (7.2), 
consider the hyperplane V perpendicular to e = (1, ..., 1) G C" 

C"Dy = {(z,e) = 0} 

and the isotropy subgroup Ge G U{N) of Je, i.e., the f/'s such that U^e = e, 
thus preserving V. That Jg = — / follows at once, from 

V^^ J l<i<n 

Since Ge stabilizes /e, there exists a unitary matrix U E Ge, diagonalizing 
X, having the property 

[UXU-\UYU-^] = [X,Y] = UleU-^ = /e, 

with 

X = diag(xi,...,x„); 

i.e., 

yXi Xj jt/ij i Oij ^ T^ Jy 

implying Y must have the form announced in (7.2). Introducing these new 
matrices into the expressions for r„ and Pnit, z) yields 

r„(t) = et^i:r*.^'det(^-x + f;H.F^-i') 

= dete^r*'^'det(-X + ti/ + 2t2l" + ...) 

n 

j=i 

andp„(t, z), as announced in (7.1) and (7.3). Note that the nrootsXj(t2,^3, ••••) 
of the characteristic equation in ti are solutions in (^2, ^3, •••) of the n-particle 
Calogero-Moser system with initial configuration coordinates {xi, ...,Xn,^i, ■■■,^n 



see T. Shiota's paper [|T5]. Thus, a solution of the discrete KP system cor- 
responds to a flag of Calogero-Moser system generated by one pair of semi- 
infinite matrices X and Y, given by (7.2), for arbitrary large n. ■ 
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Remark : Observe that for t = 0, the parameters x and z in 

PniO, z) = det{zl - Y) det (/ - {xl - Xy^ {zl - F)"^ 



are interchangeable (except for the trivial factor det{zl — Y)). This must be 
compared to the results in [|T^] and [|T^. 



8 Discrete KdV-solutions, with upper - tri- 
angular L^ 

Letting all points pi in the soliton example converge to p, all points qi converge 
to —p, and all A, converge to 1, the weights pk{z) take on the form (8.2) below. 
For future use, define the functions: 



fi = p^ sinh ^ tip^ i even 



odd 



p^ cosh ^ tip"- a odd 



odd 

and 



ge^ = p^ \z sinh ^ tip" — p cosh ^ tip" \ I even 

\ odd odd / 

= p^ \z cosh ^ tip^ — p sinh ^ tjp* ] i odd. 

\ odd odd / 



U) 



Theorem 8.1 The family of weights, 

pfc(2) = (-l)'=5W(2-p)-,5W(2 + p), forO<k<n-l, (8.2) 

leads to discrete KdV solutions, with KdV t -function^ 

n^tip' 
r„(t) = 2-e even W[foJl,...Jn-l]. (8.3) 



^Vl^[...] denotes a Wronskian with respect to the parameter p. 
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The matrix L has the property that L^ is upper-triangular, with polynomial 
eigenvectors L{t)p{t,z) = zp{t,z), given by 

W[foJu...Jn-l] 



Pn{t,z) 



•4) 



W[go,gi,...,gn~i] 

in terms of (8.1); i.e., the polynomials pn{t,z) satisfy 3-step relations of the 
following nature: 

Z^Pn{t, Z) = anPn + PnPn+1 + Pn+2- 

Then 

( oo 

>V^ = \ f = Y,a,z' such that f^''\p)- {-if f^^\-p) = {}, < A; < n - 1 



{l,z\z\...]®z{z'-pr{l,z\z\...] 



and 






Pk{u)du f 9 \ f 2p 



d 



z — u \dpj \z'^ — p"^ ^ 

k 



k = 0,...,n-l}®n^ 



1 



dp I \z^ — p2 



k = 0,...,n-l} ®n+. 



Proof : Indeed, the form of the flags Wn and W* follow from the general 
formulae (3.8) and (3.9). Therefore 

Tn{t) = det(/ijj)o<j,j<n-l 

k 



det I ( I- 
dp, 



det 



dp, 



d 



z-~^-'g{z)^^dz 
z^ — p^ 



z^-^-^g{z) 



, using (3.8) 



0<fc,jr<n-l 



dz 



z — p z + p 



0<k,j<n-l 



det I I ^ I [P- 



r(~py 



det 



y\dp. 



-J-ieE*^p'_(-p)"-i-ieS* 

0<fc,jr<n-l 
odd _|_ ( — X^^-ie odd 



." ^ gcvcn 



V 



0<A:,j<n-l 



n^tip"- 
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which is formula (8.3). 

In order to express the wave vector, one needs to compute 

\^ ' / 0<fe,jr<n-l 

= det — p"-J-ieeven eodd (^ _ p) + (_i)"-Je odd {z + p) 

nJ2Up' 
= 2"e even W[go, Ql, ..., Qn-l], 

from which formula (8.4) follows. 

Also notice that from the form of VV„, we have 



0<fc,jr<n-l 



z'Wn C Wn and thus z'W^ C Wl 
It shows that the r„(t)'s are KdV r-functions. This fact, combined with 
W* = (span{p^,p*i,...,p^_J)^ = span{p„(t,2;),p„+i(t,2;),...} C n+, 
leads to the 3-step relation: 

Z^Pnit, Z) = anPn + PnPn+1 + Pn+2, 

establishing the upper-triangular nature of L^. ■ 

Remark : Letting p ^ in p~"*^""'"^^/^r„(t) leads to the rational KdV 
solutions, i.e., the Schur polynomials with Young diagrams of type u = (n, n— 
!,...,!). 
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